Unveiling a two-dimensional electron gas with universal subbands 

at the surface of SrTiOs 

A. F. Santander-Syro 1 - 2 '*, 0. Copie 3 ' 4 , T. Kondo 5 , F. Fortuna 1 , S. Pailhes 6 , R. Weht 7 - 8 , 
X. G. Qiu 9 , F. Bertran 10 , A. Nicolaou 10 , A. Taleb-Ibrahimi 10 , P. Le Fevre 10 , G. Herranz 11 , 
M. Bibes 3 , Y. Apertet 12 , P. Lecoeur 12 , M. J. Rozenberg 13 ' 14 & A. Barthelemy 3 . 

1 CSNSM, Universite Paris-Sud and CNRS/IN2P3 Bdtiments 104 et 108, 91405 Orsay cedex, France. 

2 Laboratoire Photons Et Matiere, UPR-5 CNRS, 
ESP CI, 10 rue Vauquelin, 15231 Paris cedex 5, France. 
3 Unite Mixte de Physique CNRS/Thales, Campus de I'Ecole Poly technique, 
1 Av. A. Fresnel, 91161 Palaiseau, France and Universite Paris-Sud, 91405 Orsay, France. 
4 Universitdt Wurzburg, Experimentelle Physik VII, 91014 Wurzburg, Germany. 
5 Ames Laboratory and Department of Physics and Astronomy, Iowa State University, Ames, IA 50011. 

6 Laboratoire Leon Brillouin, CEA-CNRS, 
CEA-Saclay, 91191 Gif-sur-Yvette, France. 
7 Gerencia de Investigacion y Aplicaciones, CNEA, 1650 San Martin, Argentina. 
s Instituto Sdbato, Universidad Nacional de San Martin, CNEA, 1650 San Martin, Argentina. 
9 Institute of Physics and National Laboratory for Condensed Matter Physics, 
Chinese Academy of Sciences, Zhongguancun nansanjie 8, Beijing 100190, China. 
10 Synchrotron SOLEIL, L'Orme des Merisiers, 
Saint- Aubin-BP48, 91192 Gif-sur-Yvette, France. 
^Institut de Ciencia de Materials de Barcelona, CSIC, 
Campus de la UAB, 08193 Bellaterra, Catalunya, Spain. 
l2 Institut d'Electronique Fondamentale, Universite Paris-Sud, Bdtiment 220, 91405 Orsay, France. 
13 Laboratoire de Physique des Solides, Universite Paris-Sud, Bdtiment 510, 91405 Orsay, France. 
14 Departamento de Fisica, FCEN-UBA, Ciudad Universitaria, Pab. 1, Buenos Aires (1428), Argentina. 

* e-mail: andres.santander@csnsm.in2p3.fr 



1 



Similar to silicon that is the basis of conventional electronics, strontium 
titanate (SrTiOs) is the bedrock of the emerging field of oxide electronics 
SrTi0 3 is the preferred template to create exotic two-dimensional (2D) phases 
of electron matter at oxide interfaces {3- 5], exhibiting metal-insulator transi- 
tions 

HQ, 

superconductivity [8j, [9j, or large negative magnetoresistance [10]. 
However, the physical nature of the electronic structure underlying these 2D 
electron gases (2DEGs) remains elusive, although its determination is crucial 
to understand their remarkable properties jlll . \l2 \. Here we show, using angle- 
resolved photoemission spectroscopy (ARPES), that there is a highly metallic 
universal 2DEG at the vacuum-cleaved surface of SrTiO, 3 , independent of bulk 
carrier densities over more than seven decades, including the undoped insulating 
material. This 2DEG is confined within a region of ~ 5 unit cells with a sheet 
carrier density of ~ 0.35 electrons per a 2 (a is the cubic lattice parameter). 
We unveil a remarkable electronic structure consisting on multiple subbands 
of heavy and light electrons. The similarity of this 2DEG with those reported 
in SrTi0 3 -based heterostructures |(3, [l3] and field-effect transistors 0, I14} ] 
suggests that different forms of electron confinement at the surface of SrTi0 3 
lead to essentially the same 2DEG. Our discovery provides a model system for 
the study of the electronic structure of 2DEGs in SrTi0 3 -based devices, and a 
novel route to generate 2DEGs at surfaces of transition-metal oxides. 

In the cubic phase, stoichiometric SrTi0 3 has an empty £23 conduction manifold 
composed of three dispersive three-dimensional (3D) bands that are degenerate at the T 



point 15 1 . As schematically shown in figure 1(a), its band structure along one direction, 
say k y , consists of a weakly dispersive (heavy mass) band and a pair of degenerate strongly 
dispersive (light mass) bands. They arise, respectively, from the small and large overlaps 
along y of neighbouring Ti 3d xz , 3d xy and 3d yz orbitals, as depicted in figure 1(b) for the 
3d xy orbitals. Thus, the d xy -\ike band will be light in the xy plane and heavy along z, while 
the d yz band will be heavy along x and light in the yz plane, and the d xz band will be heavy 
along y and light in the xz plane. The resulting three Fermi surface sheets, projected on 
the xy plane of the sample surface, are shown in figure 1(c) for several bulk dopings. 

ARPES is a powerful technique to probe the electronic structure of materials |l6(. Previ- 
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ous ARPES experiments addressed the bulk electronic structure of doped SrTiC-3, revealing 
dispersing quasiparticle peaks, in-gap features, and discussing possible polaronic effects near 



the Fermi energy (Ep) 
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22|. As we shall see, our ARPES data do not conform to the 



expected lightly-doped bulk band picture, but reveal that a novel 2DEG with a complex 
subband structure, schematized in Figs. l(d,e), is realized at the low-temperature vacuum- 
cleaved surface of this material. 

Crucial to our discoveries is the comparison of samples with very different bulk carrier 
densities of n 3 p> < 10 13 cm -3 (essentially non-doped), n 3 D ~ 10 18 cm -3 (low-doped) and 
n 3 E> ~ 10 20 cm -3 (highly-doped), determined by bulk-sensitive techniques (see Methods). 
Their large doping difference is already revealed by their photographs, shown in Figs. 2(a-c): 
while the non-doped sample is transparent, the highly-doped sample is shiny black. 

Our first striking observation, shown in figure 2(d), is that even for the transparent 
non-doped sample [Fig. 2(a)], for which no bulk bands at E F are expected, the cleaved 
surface yields intense strongly dispersive bands across the Fermi level. They correspond to 
a large density of mobile carriers (see later) ensuring that, despite the electron emission, 
there is no charging of the surface of this otherwise highly insulating sample in the bulk. 
Equally stunning is the fact, shown in figures 2(d-f), that for all the studied bulk dopings the 
electronic states around Ep, and in particular the bandwidths (W = Ep— band bottom) and 
Fermi momenta (kp) of the observed bands, are essentially identical. All these observations 
contrast sharply with the large differences that would be expected from electron-doping in 
the bulk [Fig. 1(c)], where the 3D density of carriers scales as kp. The immediate conclusion 
is that the observed bands correspond to 2D states (see also the Supplementary Material). 

The spectra in figures 2(d-f), taken around r 102 and r 012 [B-like points in Fig. 2(j)], clearly 
show two strongly dispersing bands, termed 'upper and lower parabolic bands'. Through 
a closer examination, we identify other two weakly dispersing bands, termed 'upper and 
lower shallow bands'. Their presence is revealed either by rotating the light polarization 
at right angles with respect to the measurement direction, as illustrated for the non-doped 
sample in figure 2(g), or by measuring around a T-point at a larger emission angle in a 
different Brillouin zone [r^, or point C in Fig. 2(j)], as shown for the low-doped sample in 
figures 2(h-i). Note furthermore, from figures 2(d, g) that the changes in light polarization 
affect differently the two strongly dispersing parabolic bands -completely suppressing the 
upper parabolic band in the case of figure 2(g), while the increase in emission angle [compare 
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Fig. 2(e) to Figs. 2(h, i)] enhances in the same asymmetric way the negative-/;; intensities 
of the two shallow bands and weakens the parabolic bands. These dichroic effects, due to 



dipole-transition selection rules 16|, ll8| , unambiguously show that the two parabolic bands 



have different symmetries, which are also different from a third symmetry common to the two 
shallow subbands (see also the Supplementary Material). Figure 3 illustrates the complete 
set of four bands overlaid on the data. 

Let us now analyze quantitatively this band structure. The upper and lower parabolic 
bands [Figs. 2(d-f)] have bandwidths of 210 meV and 100 meV, and Fermi momenta of 
0.21 A -1 and 0.13 A 1 , respectively. Thus, they correspond to light carriers with effective 
masses m* f» 0.7m e along k y (m e is the free electron mass). The upper shallow band 
[Figs. 2(g, h)] has a bandwidth of 40 meV and a Fermi momentum kp ~ 0.3 — 0.4 A -1 , 
corresponding to heavy carriers with m* ~ 10 — 20m e . The lower shallow band [Fig. 2(j)] 
disperses from about —160 meV at T to about —120 meV at the Brillouin-zone boundary, and 
is thus fully occupied. All the above figures differ by less than 10% between the lowest- and 
highest-doping samples. As shown in figure 4, the lower parabolic band generates a circular 
Fermi surface of radius 0.21 A -1 while the upper parabolic band and the upper shallow 



band form ellipsoids along k x and k y , respectively, of semi axes 0.13 A 1 and 0.3 — 0.4 A -1 , 
yielding a 2D 4-pronged Fermi surface. 

From the area enclosed by each Fermi surface (Ap), the corresponding 2D carrier density 
is ri2D = Ap/(27i 2 ). This gives 7 x 10 13 cm -2 , or about 0.11 electrons per a 2 , for the circular 
Fermi surface, and ~ 6.2 — 8.3 x 10 13 cm -2 , or about 0.1 — 0.13 electrons per a 2 , for each 
of the elliptical Fermi surfaces. We then obtain ~ 0.31 — 0.36 electrons per a 2 (or about 
2 x 10 14 cm -2 ) coming from all the observed bands that cross Ep. 

We now rationalize the observed electronic states as the subbands of a 2DEG confined 
within a few unit cells at the surface of SrTiCV To this end, we consider a potential 
Vo at the surface that confines the electronic motion along z [inset of Fig. 1(d)]. This 
lowers the energy of the bands by about Vo, similar to the 'band-bending' at semiconductor 
heterostructures, and produces an energy splitting between the different eigenstates that 
is inversely proportional to their effective masses along z (m*). The resulting subband 
structure, depicted in Fig. 1(d), consists on a single d xy -\ike band and two d xz /d yz bands 
that are degenerate at T. As the d xy band has a very large m*, the attractive confining 
potential will merely pull it beneath Ep (its energy-splitted eigenstates will have a negligible 
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separation). Thus, we identify this band with the lower parabolic band in our spectra of 
figure 2, and note it E\(d xy ). On the other hand, the d xz and d yz subbands, light along z, 
will show large energy splittings. They are noted E n {d xz / yz ) (n = 1, 2, . . . ) in figure 1(d). 



Other effects beyond 
temperature tetragonal 



;his simplified analysis, such as spin-orbit coupling and/or the low- 
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24| and possibly orthorhombic 23] distortions, can lift the de- 



generacy between the d xz and d yz subbands, as represented in Fig.l (e). This would induce 
a coupling, resulting in hybridization between the light and heavy subbands, as indeed 
evidenced by our data for the lower parabolic band and the shallow bands (Fig. 3, LV 
polarization). Another possibility would be a surface reconstruction, although this is not 
evidenced by our data -which follows the periodicity of the unreconstructed bulk lattice 
without apparent band folding. Resolving the detailed origin of such degeneracy lift, which 
is not relevant to the discoveries reported here, will require future investigations. Hence, we 
identify the lower shallow band and upper parabolic band in our spectra (Fig. 2) with the 
splitted doublet Ei(d xz ) and Ei(d yz ), respectively. From the data, the doublet splitting is 
A « 60 meV. We finally identify the upper shallow band with E2(d xz ). Given that the bot- 
tom of this band is at —40 meV, its upper partner E2(d yz ) would occur at about +20 meV, 
and cannot be observed. 

This subband hierarchy respects the symmetry considerations previously deduced from 
the effects of light polarization and emission angle, rendering additional support to the 
picture of a confined 2DEG. 

To characterize quantitatively this 2DEG, we use a simple schematic model. We assume 
that the surface has a homogeneous positive charge density, generating a confining electric 
field F inside the solid, and thus a triangular wedge potential of depth Vo at the surface 
(z = 0) and value V(z) = — Vo+eFz along z, as depicted in the inset of Fig. 1(d). The details 
of the calculations are presented in the Supplementary Information. From the experimental 
splitting E2(d xz ) — E\{d xz ) = 120 meV between the first two subbands of the d xz orbitals, 
and their light mass m* w 0.7m e along z, we deduce a strong confining field F w 83 MV/m. 
From Ei(d xy ), the confining potential is estimated to be Vo ~ —260 meV. The deduced 
values of Vo and F, and the measured effective masses and doublet splitting of 60 meV 
between the d xz and d yz subbands, suffice to retrieve all the observed subband energies. The 
width L of the confined 2DEG results from the spatial extension of the highest occupied 
state E 2 (d xz ), or independently from the approximate expression eFL E 2 (d xz ) — Ei(d xz ) 
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[inset of Fig. 1(d)], yielding consistent values in the range L 14.5 — 18 A rj 4 — 5 unit cells. 
Additionally, from the confining field F and the polarizability of the medium, we calculate 
an induced surface charge density of ~ 0.25 electrons per a 2 , in good agreement with the 
experimental value. 

The band-bending due to confinement should induce an energy down-shift of about Vq 
of the 0-2p valence band, creating a gap with respect to Ep that is larger than the optical 
gap. This effect has been reported in previous experiments |18l. |20| . and is also observed in 
our data (Supplementary Material). Aiura and coworkers [18] demonstrated that it is due 
to surface oxygen- vacancies, which appear thus as the most reasonable candidate to explain 
the origin of the 2DEG observed in our experiments. In fact, cleaving is likely to produce 
a massive removal of surface oxygen, much larger than and basically independent of the 
bulk concentration of vacancies, providing 2 dopant electrons per created vacancy. These 
two electrons will delocalize within the potential wedge created by the positively charged 
layer of surface vacancies. Note that cleaving this almost cubic oxide system necessitates 
overcoming the strong binding electric forces between atoms, equivalent to inducing a 'me- 
chanical dielectric breakdown'. Interestingly, from our data analysis, the large electric field 
generated at the surface of SrTiC>3 is of the order of typical dielectric breakdown fields 25] , 
and thus compatible with a cleaved surface in this material. 

To give additional support to our picture and data analyses, we performed ab-initio 
calculations including oxygen vacancies at the surface. The results are consistent with the 
above simplified wedge-potential model (Supplementary Material). 

Strikingly, the sheet carrier density, the confinement size, the confining field F and po- 
tential Vo, and the presence of both mobile and quasi-localized carriers in this novel electron 
gas at the surface of SrTiC>3 compare well with the characteristics of other 2DEG's at dif- 
ferent types of SrTi0 3 interfaces (], [§], [9], Ij], 13, 26-28]. All these remarkable consistencies 



suggest that all these 2DEG's may be understood on a common basis. In particular, note 
that different electron confinement mechanisms, such as the electric field created by a polar 
LaA10"3 layer j^} or the direct application of an electric field 9], would recreate a 2DEG 
confined within a nanometric layer of SrTiOs with a subband electronic structure similar in 
nature to the one directly unveiled by our measurements. 

Our experiments directly show that the Ti-3d xy states of light and highly mobile carriers 
become the first available levels for the 2DEG at the surface of SrTiOs. Altogether, our 
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results demonstrate the fundamental paradigm for understanding the properties of oxide 



surfaces and heterostructures, namely the 'electronic reconstruction' into several partially- 
occupied subbands at the near-interface region 
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29 



30j. Our discoveries also raise 



the appealing possibility that oxygen vacancies may be induced to aggregate or self-organize 
at the surface of SrTi03 or related interfaces -a situation that might be at the origin of 
the reversible patterning of conducting nano-lines in LaAlOs/SrTiOs interfaces 

aa 

. In a 

similar way, nano-patterning of oxygen vacancies at the surface of SrTiOs could provide a 
model system for low dimensional confined electron gases. More generally, we suggest that 
the engineering of atom- vacancies at the surfaces of transition metal oxides will open a new 
avenue in correlated-electron surface science. 



METHODS 

Sample preparation and measurement technique 

High-quality single crystals of SrTiOa were doped with oxygen vacancies by high- 
temperature treatment in low-oxygen pressure, except for the non-doped sample, which 
was a bare SrTi03 substrate measured without any previous preparation. The bulk carrier 
concentrations were deduce from Hall and resistivity measurements. 

The ARPES experiments were done at the Synchrotron Radiation Center (SRC, Univer- 
sity of Wisconsin, Madison) and at the Synchrotron SOLEIL (France), using 47 eV (SRC) 
and 45 eV (SOLEIL) linearly polarized photons, a Scienta 2002 detector with horizontal slits 
for the highly-doped sample (SRC), and a Scienta R4000 with vertical slits for the low-doped 
(SRC) and non-doped (SOLEIL) samples. The above photon energies, close to resonant pho- 



toemission at the Ti 3p — > 3d edge, yield intense quasi-particle peaks in SrTiOs [18|. The 
momentum and energy resolutions were 0.25° and 20 meV, respectively. The diameter of 
the incident photon beam was smaller than 100 /im. The samples were cleaved in-situ along 
the c-axis at 20 K (SRC) and 10 K (Soleil), in pressure lower than 6 x 10~ n Torr. For each 
set of measurements, one narrow highly-emitting terrace was kept fixed. The results have 
been reproduced in at least two cleaves for each bulk doping. 
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Measurement geometries 

The sample's surface defines the xy-pleme. The measurements on the /c-space hemisphere 
were taken along (010) (or k y ). The data in Figs. 2(d, g) and Fig. 3 were collected in B- 
like points slightly above r 10 2, with the photon polarization along y (LV) in Fig. 2(d), and 
along the rrz-plane (LH) in Fig. 2(h). The data in Figs. 2(e, f) [Figs. 2(h, i)] correspond to 
B-like (C-like) points slightly above r 012 (slightly below r 112 ). For these points, the photon 
polarization is not parallel to any of the sample's symmetry directions or planes, and the 
spectra have different symmetry mixtures. The data in Fig. 4 were collected with LH photon 
polarization. The Supplementary Information discusses further the measurement geometries 
and their corresponding dipole selection rules. 
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FIG. 1: Schematic electronic structure of SrTi03 and effects of electron confinement, a, Bulk 

conduction band of SrTi03 along k y , consisting on a heavy d xz band (green) and a doublet of light d xy — d yz 
bands (red and blue). These bands stem from the small and large overlaps of the Ti-3d orbitals along the 
y direction, as depicted in b for the case of 3d xy orbitals. c, Resulting 3D Fermi surfaces (cut along the 
rry-plane) for several bulk dopings. The last panel shows, for comparison, the 2D Fermi surface for a 2DEG 
of density ~ 10 cm" 2 (hybridization between different Fermi-surface sheets is not included), deduced from 
a tight-binding model [12 1. Colors indicate the character of each Fermi surface sheet along k y . d, Quantum- 
well states, or subbands, resulting from the confinement of electrons near the surface of SrTiC>3. The inset 
shows a wedge-like potential created by an electric field of strength F at the surface, used as a simple model 
to analyse the ARPES data (main text). The size L of the confined 2DEG can be estimated from the 
extension of the highest occupied state (n — 2 in our case) or approximately from the energy difference 
between the lowest (n = 1) and highest occupied states, e, Additional degeneracy lifts at T occur due to 
spin-orbit coupling, tetragonal and orthogonal distortions, or possible surface reconstructions. 



12 





FIG. 2: Universal electronic structure at the surface of SrTiC>3. a-c, Photographs of the samples 
studied by ARPES. Their bulk dopings are written beneath, d-f, Corresponding energy-momentum inten- 
sity maps close to the T point (B-like points in panel j). g, Effects of changing the photon polarization 
(specifications below) from linear vertical (LV, panel d) to linear horizontal (LH), illustrated here for the 
non-doped sample. The upper shallow band appears, and the upper parabolic band disappears. The effects 
of collecting electrons around a T-point at a larger emission angle in a different Brillouin zone (point C in 
panel j) are also shown, for the case of the low-doped sample, in the energy- momentum intensity map of 
panel h, and on its second derivative (only negative values shown) in panel i, where the dispersion of the 
lower shallow band is indicated by a dashed red line. Here, the negative-fc intensities of the shallow bands 
are enhanced, while the intensities of the parabolic bands decrease, j, 2D representation of the 3D Brillouin 
zones of SrTiOs (dashed lines) in relation to the hemispherical measurement surface in momentum-space at 
hv = 47 eV lallSl- The measurement temperature (10 K or 20 K) is indicated in each panel. 
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10K LV 



10K LH 




FIG. 3: Summary of subbands for the 2DEG at the surface of SrTiC>3. Side-by-side representation 
of the LV and LH spectra around Tio2 for the non-doped sample. An elongated aspect ratio has been used 
to better discern the shallow bands. The dotted lines are tight-binding representations of the bands 12[, 



following the same color scheme as in Figs. l(d, e). Hybridizations between the lower parabolic band and 
the shallow bands are observed on the spectra taken with LV polarization. 
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FIG. 4: Fermi surface of 2DEG at the surface of SrTiC>3. a, ARPES Fermi surface map (integrated 
within ±20 meV around Ep ) across several Brillouin zones for a photon energy hv — 47 eV, displayed 
on the spherical shell of 3D fc-space and projected on the k x — k y plane. The Fermi-surface circle from 
the lower (d^y-like) light parabolic band is best observed in the Foo2 and Foi2 (and equivalent) Brillouin 
zones, while the two elliptical perpendicular Fermi surfaces from the upper parabolic band (d yz ) and the 
upper shallow band (d xz ) are best observed off-normal emission (Foi2, Fn2 and equivalent), b, Second 
derivative (negative values only) of the Fermi surface map over Fio2, showing portions of the circular and 
the horizontal elliptical Fermi surfaces. The intensity of the vertical ellipse is very weak here due to dipole 
selection rules, c, Collected Fermi momenta from high-statistics measurements around Tio2 (white circles) 
and Tii2 (grey circles), and comparison to the Fermi surfaces resulting from the tight-binding bands of 
figure 3 (colors indicate the character along k y and follow the scheme of precedent figures). All panels 
correspond to measurements on the low-doped sample. 
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Abstract 

The angle-resolved photoemission spectroscopy (ARPES) data on SrTiOs show the existence 
of a two-dimensional electron gas (2DEG) beneath the surface of this material. We provide here 
additional evidence of its 2D-like nature, showing that the bands do not disperse along the k z [or 
(001)] direction, although their intensity is strongly modulated as a function of photon energy, 
photon polarization and emission angle for different Brillouin zones. We also provide the details 
of the modeling of these bands in terms of a confining wedge-like potential well. Additionally, we 
show that in our data the 0-2p valence band is shifted to larger binding energies with respect 
to what is expected from the bulk band gap ('band bending'), by an amount comparable to the 
confining potential depth. We finally present ab-initio calculations that support the picture of the 
formation of a 2DEG due to oxygen vacancies at the surface of SrTiO-3. 
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Photon energy dependence of the electronic subbands at the surface of SrTi03 



Figure SF [T] presents data obtained at normal emission as a function of photon energy. 
The data were collected using linear horizontal (LH) light polarization (polarization vector 
in the xz plane), and detection along x [horizontal (H) slits] in the sample's xz-plane, which 
in this geometry also coincides with the light incidence plane. The final electron state has 
to be even with respect to the sample's mirror plane to be detected by the analyzer (the xz- 
plane in this case) -otherwise the electron wavefunction would be zero in this plane. Hence, 
the dipole-transition selection rules imply that the initial state and the dipole operator 
(equivalently, the polar ization vector) must have the same symmetry with respect to the 



mirror plane 



SR3 



SR4J |. In the geometry of figure SF[H the light polarization is even under 



reflection on the xz plane. Thus, only initial states of even parity with respect to this plane 
can be detected. Hence, only the upper d^-like band, which is furthermore light along x, is 
observed in this geometry (experimentally, the lower d xz -like band has a very weak intensity 
in this geometry). 

As observed in Fig. SF [He, d), the upper d^-like band does not disperse as a function 
of photon energy, or equivalently as a function of k z , confirming its 2D surface-like charac- 
ter. However, it presents a marked intensity modulation as a function of photon energy, a 
behav iour also o bserved in surface states and quantum well states of simple metallic sys- 
tems 



SRI 



SR2| . Thus, for instance, while the band is clearly seen at hu = 30 eV, its 



intensity is very weak around hu = 47 eV, where the rest of our experiments were per- 
formed. As explained in the main text, and also discussed below, at this photon energy the 
d xz -like bands and the other subbands are distinctly observed in other Brillouin zones and 
other polarization/slits geometries. 



Intensity modulations of the electronic subbands with emission angle and photon 
polarization 

Fig. SF|2]shows the energy-momentum intensity maps for cuts at three points in reciprocal 
space labeled as A, B and C. The spectra at points B and C were already discussed in the 
main text, and are again shown here for clarity and completeness in the discussion. 

The three points A, B and C correspond to different positions along the (001) (or k z ) 
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Figure SF 1: Photon-energy dependence of the ARPES spectra of SrTiOa. a, Schematics of 
the experimental geometry for the data in this figure. The sample is the dark parallelepiped, the xz mirror 
plane is transparent, the light beam is represented by the blue arrow, its polarization by the red double 
arrow, and the collected electron beam by the dashed orange arrows. The same conventions are used in the 
rest of the figures in this Supplement, b, Representation of the 3D reciprocal space. The region covered by 
varying the photon energy from 25 eV to 100 eV is shown red dashed arrow and the thick red lines. The 
open circles are T points, the dashed lines are Brillouin-zone edges, and the violet circle is the spherical shell 
on which the measurements at hv = 47 eV are carried out. To convert the photon ene rgy d ependence into 



SR4| . c, Integrated 



k z dependence, an 'inner potential' of 12 eV was used, as deduced in previous works 
photoemission spectra at norm al emission over a wide energy range, showing the oxygen-vacancy related 
intra-gap states at 1.3 eV 



SR4 ] and the quasi-particle (QP) peak near Ep. Notice the large variations of the 
QP peak intensity as a function of photon energy (varying from 25 eV to 100 eV in steps of 5 eV). d, Zoom 
over the QP peak in c, after normalizing by the area of the peak. The peak does not disperse with photon 
energy (26 — 100 eV in steps of 2 eV), indicating its 2D character, e, Examples of two energy-momentum 
intensity maps at hv = 30 eV, at roo2> and hv = 46 eV, close to the zone edge. The band corresponds to 
the d xz states. Its intensity is almost suppressed at hv — 46 eV, but the dispersion is unchanged. 
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axis and different emission angles, as represented in Fig. SF [2(a) [SR4j . |SR6[ |. such that 
r 002 A = 0.7(f) = 0.56 A -1 , r i 2 B = 0.25(f) = 0.2 A" 1 and r 112 C « -r i 2 B (a is the cubic 
SrTi03 lattice parameter). All the spectra in this figure were measured at hv = 47 eV with 
LH light polarization and detection along y [vertical (V) slits], as schematically shown in 
Fig. SF |2(b). At normal emission (point A), the detection plane coincides with the sample's 
yz-p\ane. Exactly at T, as explained in the previous section, the only symmetry-allowed 
states are d xz -like, which have however a very weak intensity at this photon energy. Note 
that in this case the polarization vector has components both along the x and z directions, 
the first being odd with respect to the detection plane, and the second even with respect 
the same plane. Thus, away from T, at emission along the yz-plane, both d^-like bands 
(odd with respect the yz-plane) and o^-like bands (even with respect the same plane) can 
be observed. 

On the other hand, for points B-like and C-like, the photon polarization is not parallel to 
any of the sample's symmetry directions or planes, and the spectra have different symmetry 
mixtures depending on the measurement point in fc-space (see also further). 

From figure SF [2(c), we see then that at normal emission (point A) the two strongly 
dispersing parabolic bands are distinctly observed. These bands remain unchanged at points 
B and C [Figs SF |2(d,e)], though their spectral intensity decreases as the emission angle 
increases. Note, in particular, that although points B and C are equidistant from T, the light 
parabolic bands are more intense at point B, where the emission angle is lower. Conversely, 
the two weakly dispersing bands (the 'upper and lower shallow bands') are best seen at 
point C, where the emission angle is largest. Figures SF[2(f, g) show the energy distribution 
curves (EDCs) at point C integrated over intervals of 0.1 A -1 . The lower shallow band is 
clearly seen in the raw data dispersing from V (blue EDC) to the zone edge (green EDC). 
Additionally, as noted in the main paper, the same asymmetric enhancement of the negative- 
k intensities for the two shallow bands is observed in point C. This indicates that the two 
shallow bands have the same symmetry. These shallow bands rapidly loose intensity as the 
emission angle decreases, to completely disappear in the spectra at normal emission, further 
indicating that their symmetry (or more generally, their spectral function) is not the same 
of either of the parabolic bands. 

From the data of Fig. SF [21 it is then apparent that there is a unique underlying band 
structure at all three points A, B and C, further proving the 2D-like character of all the 
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Figure SF 2: Intensity modulations of the electronic subbands as a function of emission angle. 

a, Schematic representation of the reciprocal space of SrTi03 over several Brillouin zones in the (100) — (001) 
and (110) — (001) planes. The continuous line is the spherical shell on which the measurements at hv — 47 eV 
are performed, b, Schematics of the experimental geometries for the data taken at points A, B and C in 
this figure, c-e, Energy-momentum intensity maps near Ep at points A, B and C represented in panel a. 
f-g, EDCs at point C integrated over intervals of 0.1 A -1 . The blue EDC is located at T, and the green 
EDC is at the zone edge 

observed bands. 

As described in the main text, the photoemission intensity from the different subbands 
is also strongly modulated by the photon polarization. For the sake of completeness in the 
discussion, we show again in Fig. SF|3]the energy-momentum intensity maps and EDCs taken 
at Tio2 (point B) for the non-doped sample. These data were collected at hv = 45 eV with V- 
slits (detection along y) using both LV and LH polarized photons, as shown in Figs. SF [3^a, 
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b). In this geometry, where electrons are detected off-normal emission, only the xz-plane 
(coinciding with the light-incidence plane), where Tio2 lies, is a high-symmetry plane. For LV 
polarization, the light electric field is odd with respect to this mirror plane. Thus, only the 
d xy - and d yz -Y\ke bands, which are both odd with respect to the xz-plane, can be detected, 
as indeed observed [Figs. SF[3^c, d)]. On the other hand, for LH polarized light, the electric 
field is even with respect to the mirror plane, and only d xz -like bands can be detected at 
T102, also in agreement with the observation of the upper shallow band [Figs. SF [3te, f)]. 
Away from Fio2 along the detection direction there is symmetry mixing. Experimentally, 
the lower parabolic band, but not the upper parabolic band, is again observed. 



Band bending of the 0-2p valence band 



Figure SFH^a) shows the angle-resolved data over a wide energy range at normal emission, 
LH polarization, V-slits and hv = 47 eV [point A in figure SF EJ^a)] for the sample that 
has bulk doping n 3 £, = 10 18 cm -3 . The dispersion of the oxygen valence band, roughly 
located in between -9 eV and -4 eV, is clearly observed, attesting the high-quality of the 
surfaces obtained after cleaving. As seen from figures SF HJ^a, b), the peak position of 
the valence-band maximum (VBM) is located at E ~ —4.35 eV, and the leading-edge of 
the valence band is located at about —3.9 eV. These values are systematically lower by 
200-500 meV than the pe ak po sition of the bulk VBM, predicted to be about —3.9 eV by 
tight-binding calculations |SR6j |. and the bulk optical band gap of 3.75 eV [SR7| . shown for 
compa rison . Thes e results are consistent with the dispersions and shifts reported in previous 
works 



SR4 



SR6| . More important, the observed energy down-shift of the oxygen valence 



band is in quantitative agreement with the confining potential at the surface, estimated to 
be Vq ~ —300 meV (see below and main paper), and responsible for the formation of the 
two-dimensional electron gas and the subbands beneath the surface of cleaved SrTiCV This 
indicates that such a band shift of the 0-2p valence band is d ue to band-bending by the 
above confining potential. Furthermore, as noted in Ref. SR4J, when the cleaved surface 
is exposed to oxygen, the 0-2p valence band shifts back by about 200 meV towards Ep, 
lending support to the hypothesis that the confining potential, in the case of cleaved SrTiC>3, 
is due to surface oxygen vacancies. 
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Figure SF 3: Intensity modulations of the electronic subbands as a function of photon polar- 
ization, a, b, Schematics of the experimental geometries for the data taken at point B [Fig. SF^a)] in this 
figure, c, d, Energy-momentum intensity map and EDCs for LV polarized light, e, f, Same measurements 
for LH polarized light. Blue EDCs are located at Tio2- Both sets of spectra were obtained in the same 
spatial spot of the sample, at T = 10 K. 
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Figure SF 4: ARPES spectra and band-bending of the 0-2p valence band in SrTiOa. a, Second 
derivative of the energy-momentum intensity map over a wide energy range, showing the dispersion of the 
0-2p valence band at energies between 4 — 9 eV. The non dispersive oxygen- vacancy- induced intra-gap states 
at —1.3 eV, and the QP band ar Ep, are also observed, b, Raw EDC at normal emission (k v = 0) from the 



SR7I ] is shown 



data in a (blue curve). The energy-dependent optical extinction coefficient measured in Rcf. 
for comparison (black curve). The leading edge of the photoemission spectrum is shifted down by about 
300 meV with respect to the optical absorption leading edge. 



Wedge model 



We approximate the effect of the confining potential acting in the ^-direction by a poten- 
tial wedge V(z) = V^ + eFz, where e is the charge of the electron and F is the strength of the 
ele ctric field a long z. The quantized eigenenergies E n are given in very good approximation 



by 



SR8. 



SR9|: 



V + 



h 2 



2m* 



1/3 r 



3tt 

T 



l 



n 



eF 



2/3 



(1) 



where n = 1,2,3,. . ., and m* is the effective mass in the ^-direction. The exact eigenenergies 
have n — \ in this equ ation replaced by 0.7587, 1.7540, and 2.7575, respectively, for the 
three lowest solutions [SR9J. The average value of z for the n-th subband is eF (z n ) = 



l\E n -V \ [SR9|. 



Replacing with the appropriate physical constants, one gets 

= (n-i) F 2/3 , (2) 

where m e is the free electron mass, E n and Vb are expressed in eV, and F is expressed 
in V/m. This expression allows to determine the field strength F from the experimental 
subband splitting AE = E2 — E% between the first to subbands of d xz / yz character. Notice 
that, in cubic symmetry, the d xz - and d yz -\Ske. bands have the same small effective mass in 
the ^-direction, with m* pa 0.7m e from our experiments. 

From the experimental value AE pa 120 meV, we get F pa 83 MV/m. 

Given the wedge-like form of the potential, the states with higher energy will have the 
largest extension along z. In our case, the observed highest occupied state is the upper 
shallow subband, which would correspond to the n = 2 state of the d xz -like band. Thus, 
with the above field strength, we can estimate the depth (L) of the confining potential from 
the approximation eFL pa [E 2 {d xz ) — Ei(d xz ] (note that Vb is not needed to compute this). 
This gives L pa 14.5 A, or about 3.7 unit cells. 

The potential at the bottom of the wedge well (Vb) can, in turn, be estimated from the 
field strength and the bottom of the d xy parabolic band [Ei{d xy ) pa —210 meV], which has 
a heavy mass along z of m* pa 20m e . This yields Vq pa —260 meV. This value allows an 
independent estimate of the width of the well from the average value of z for the n = 2 
subband, namely eF (z 2 ) = \\E 2 — V \. Using E 2 {d xz ) = —40 meV, we obtain (z 2 ) = 17.7 A, 
or about 4.5 unit cells. Additionally, as the potential V{z) remains negative {i.e., confining) 
while eFz < Vb, this yields L max pa 31 A pa 8 unit cells as a maximum bound for the width 
of the 2DEG. These independent estimates provide important crosschecks for the internal 
consistency of our analysis. 

The energy of the bottom of the different d yz / xz subbands is then readily calculated from 
their light effective masses along z and the values of Vb and F. One obtains: Ei{d yz / xz ) = 
— 100 meV, and E 2 {d yz / xz ) = +20 meV. This agrees with the observations: only the first 
d yz subband (the small light parabola in the xy-plane) lies below Ep. The experimentally 
observed degeneracy lift between the d xz and d yz bands is not taken into account by this 
simple model. From the data, we know that the Ei{d xz ) subband (the lower shallow subband) 
is shifted by about 60 meV below its d yz partner. One then expects that the upper shallow 
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— 60meV = —40 meV, again in excellent 



subband lies at an energy E2{d xz ) = E 2 {d xz / yZi 
agreement with the observations. 

Finally, we can also use our estimate of F to compute the superficial charge density a 2D 
induced by the confi ning; potential. We follow the same type of calculation used by Ueno 
and coworkers (Ref. SR8[, Supplementary Material), namely: 

-F 



2 a2D 



e e(F')dF', 



(3) 



where eo is the vacuum dielectric constant and e(F) is the material's dielectric permittivity. 
Empirically, it is well established that the dielect ric per mittivity of SrTiOs depends on the 
electric field strength F according to the relation [SRlOj: e{F) « 1/(4 x 10~ 5 + 5 x 1CT 10 F), 
where F is in V/m. We then readily obtain a ~ 0.25 e/a 2 , where a = 3.905 A is the 
cubic lattice parameter. This value is consistent with the experimental electron count ~ 
0.31 — 0.36 e/a 2 , obtained from the area of the Fermi surfaces (see main text). 



LDA slab calculations of the electronic structure in the presence of surface oxygen 
vacancies 

To illustrate the effects on the band structure of SrTiOs due to the presence of oxygen 
vacancies near the surface, we performed ab initio density functional calculations within the 
local density approximation. We adopted the Wien2k code and modelled the surface by 
considering slabs of 1 x 1 x 12 cells. We considered both SrO and Ti02 terminations, and 
put one oxygen vacancy defect at the surface. Our goal here is to understand, through a 
realistic yet computationally simple calculation, whether the presence of such vacancies can 
yield a 2DEG at the surface with the concomitant energy lowering for the slab states. Thus, 
in this simple initial study, we considered the undistorted cubic cell and did not relax the 
lattice. Our results are presented in Fig. SF [51 They show that the main features of the 
2DEG at the surface of vacuum-cleaved SrTiOs, discussed in the main text, are borne out by 
this calculation. Namely, that the bottom of the t 2g manifold around the T-point is lowered 
in energy by the effective attractive potential by as much as 200 meV. Moreover, the results 
also correctly capture the fact that the d xy parabolic band, due to its large effective mass 
along the ^-direction, is pulled down more that the d xz / yz doublet. 

While it is appealing that these generic features are well captured by the LDA results, 



11 



one should bear in mind that due to the technical limitations in the cell size that can be 
considered, our calculations remain on the qualitative level. A larger computational effort, 
including clusters of surface vacancies, is in progress. 
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Figure SF 5: Band structure calculation of 1 x 1 x 12 slabs of cubic SrTiC>3 with one oxygen 
vacancy on the surface, a, b Calculations for SrO and Ti02 terminated surfaces, respectively. A subband 
structure of levels is created due to the z-direction confinement, with a large level spacing for the d xz / yz 
doublets, which share the same light mass m*. In contrast the d xy levels are hardly split due to their large 
m*. Panels c and d include one oxygen vacancy defect at the SrO and TiC>2 surfaces. One clearly observes 
the double effect of doping and confinement. The bottom of the levels is pulled beneath the Fermi energy. 
Panels e and f show an enlarged detail of the energy region near Ep. Colors indicate the character of each 
band along k y , following the same convention as in the main paper: d xy -\ike bands are red, d yz -like bands 
are blue, and d xz -like bands are green. 
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